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1. Introduction 

Current terrestrial neutrino oseillation experiments searehing for fundamental information in 
the neutrino sector, such as the neutrino mass hierarchy and whether CP violation occurs for lep¬ 
tons, usually employ experimental designs which rely on the partial oscillation of a beam of muon 
neutrinos into electron neutrinos.[l, 2] These experiments build large detectors of heavy materials 
to maximize the rate of neutrino interactions, and then examine the energy distribution of the neu¬ 
trinos that do interact with the detector, comparing the observed spectrum with predictions based 
on hypotheses of no oscillation or oscillation with given parameters. 

Correct prediction of the observed energy spectrum for electron neutrino interactions—on 
which these oscillation results depend—^requires an accurate model of the rates and outgoing parti¬ 
cle kinematics. This, in essence, boils down to a need for precise Vg cross-sections on the detector 
materials in use. And yet, because of the difficulties associated with producing few-GeV electron 
neutrino beams, even when including very recent results, only two such cross-section measure¬ 
ments exist[3, 4]. Furthermore, the small statistics and inclusive nature of both of these measure¬ 
ments make their use as model discriminators challenging. Instead, most simulations begin from 
the wealth of high-precision cross-section data available for muon neutrinos and apply corrections 
such as those discussed in ref. [5] to obtain a prediction for Vg. 

We offer here a preliminary result in an effort to produce a higher-statistics cross-section for a 
quasi-elastic-like electron neutrino process, which is among the dominant reaction mechanisms at 
most energies of interest to oscillation experiments. We use the MINERvA detector, which consists 
of a central sampling scintillator region, built from strips of lluoror-doped scintillator glued into 
sheets, then stacked transverse to the beam axis; both barrel-style and downstream longitudinal 
electromagnetic and hadronic sampling calorimeters; and a collection of upstream passive targets 
of lead, iron, graphite, water, and liquid helium. The detector design and performance are discussed 
in full detail elsewhere. [6] MINERvA occupies space in the NuMI v^ beam, where it was exposed 
to a flux of ~ 99% v^ and ~ 1% v^ mostly between 3-5 GeV for this dataset. We also compare the 
result for Vg to a similar, previous MINERvA result for v^ to evaluate how similar they are. 


2. Signal definition 

In traditional charged-current quasi-elastic neutrino scattering, CCQE, the neutrino is con¬ 
verted to a charged lepton via exchange of a W boson with a nucleon, resulting in the following 
reaction: V/u — p. (Antineutrino scattering reverses the lepton number and isospin: V/p —> /+«.) 
Because the MINERvA detector is not magnetized, we cannot differentiate between electrons and 
positrons on an event-by-event basis. Moreover, hadrons exiting the nucleus after the interaction 
can re-interact and change identity or eject other hadrons[7]; furthermore, pairs of nucleons corre¬ 
lated within the initial state may cause multiple nucleons to be ejected by a single interaction [8, 9]. 
Therefore, we define the signal process “phenomenologically,” by its final-state particles: we search 
for events with either an electron or positron, no other leptons or photons, any number of nucle¬ 
ons, and no other hadrons. We call this type of event “CCQE-like.” We also demand that events 
originate from a 5.57-ton volume fiducial volume in the central scintillator region of MINERvA. 
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3. Event selection and backgrounds 

Candidate events are selected from the data based on three major criteria. First, a candidate 
must contain a reconstructed cone object of angle 7.5°, originating in the fiducial volume, which 
is identified as a candidafe elecfromagnefic cascade by a mulfivariafe PID algorifhm. The laffer 
combines defails of fhe energy deposifion paffern bofh longifudinally (mean dE / dx, fraction of 
energy af downsfream end of cone) and fransverse fo fhe axis of fhe cone (mean shower widfh) 
using a ^-nearesf-neighbors (kNN) algorifhm. Secondly, we separafe elections and positions from 
phofons by cuffing evenfs in which fhe energy deposifion af fhe upstieam end of fhe cone is consis- 
fenf wifh fwo particles rafher fhan one (since phofons fypically inferacf in MINERvA by producing 
an election-position pair). Af fhis poinf, fhe cone objecf becomes fhe election candidafe. Our fi¬ 
nal criferion is an affempf fo selecf CCQE-like inferacfions using a classifier we call “extia energy 
fracfion,” *P, which, when an evenf’s visible energy nof inside fhe election candidafe or a sphere of 
radius 30cm cenfered around fhe cone verfex is denofed “extia energy,” is defined as: 

(3.1) 

^electron 

Our cuf is a funcfion of fhe fofal visible energy of fhe evenf. The cuf af fhe mosf probable fofal visi¬ 
ble energy, Eyis = 0.4 GeV, is illusfrafed in fig. 1. Einally, we refain only evenfs wifh reconsfrucfed 
election energy in fhe range 1 GeV < Eg < 10 GeV. Here fhe lower bound excludes a region where 
fhe background esfimafe is sfill under furfher sfudy, and fhe upper bound resfricfs fhe sample fo 
evenfs where fhe uncerfainfies on flux predicfion are folerable. The disfribufion of evenfs selecfed 
by this sequence is shown in fig. 2. 



Extra energy ratio ^ 

Eigure 1: Sample cuf on 'P (defined in fhe fexf) af fhe mosf probable evenf visible energy, Ey/i = 0.4 
GeV. 
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Figure 2: Event sample after all selection cuts. 


As fig. 2 shows, even after the final selecfion, a significanf fracfion of fhe sample is predicfed 
fo be from background processes. We affempf fo consfrain fhe background model by examining 
sidebands in fwo of fhe variables already menfioned. The firsf of fhese is in fhe dE/dx measured 
af fhe fronf of fhe elecfron candidafe; by choosing a sample af larger values, we can obfain a side¬ 
band rich in photon background evenfs. The second sideband is in fhe exfra energy fraction 'F; 
a sample of evenfs af larger 'F consfifufes a sideband rich in inelasfic background. We use fhese 
sidebands fogefher fo fif fhe normalizations of fhe fhree major backgrounds: Vg non CCQE-like, 
non-Ve coherenf pion, and ofher inelasfic evenfs. The fhree background classes’ normalizations are 
filled simullaneously, using dislribulions in bolh reconslrucled candidafe elecfron angle and energy, 
across fhe fwo sidebands, fo obfain a single scale factor lhaf represenls fhe besl eslimale of fhe folal 
normalization of fhe background as compared to fhe prediclion from GENIE. We obfain a scale 
factor of 0.69; Ibis overall reduction is a similar Irend to lhaf observed when similar procedures 
were performed on ofher MINERvA analyses. Subsequenl to fhe conslrainl, we scale fhe back¬ 
grounds in fhe signal region and sublracl Ihem from fhe dala. We Ihen compare if fo fhe simulaled 
prediclion of fhe signal process. 


4. Cross-section result 

We calculale fhree differenlial cross-seclions in elecfron angle, election energy, and four- 
momenlum Iransfered from neutiino fo nucleus Q^, as well as fhe folal cross-section vs. neu- 
Irino energy. Eor neulrino energy and Q^, we employ fhe commonly-used CCQE approximalions 
(assuming a slalionary fargel nucleon) which allow us fo compule Ihem from jusf fhe lepton kine- 
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matics: 


^QE _ ml-inip-EbY-rnl + 2{rnp-Ei,Ee) 
^(jttp Ei^ Eg pg cos Qg^ 

Qqe = 2E^^ {Eg - pg cos ee)-m] 


(4.1) 


(4.2) 


Differential cross-sections are calculated in bins i according to the following rule for sample vari¬ 
able ^, with e representing signal acceptance, <I> the flux integrated over the energy range of the 
measurement, the number of targets (CH molecules) in the fiducial region. A, the width of bin i, 
and Uij a matrix correcting for detector smearing in the variable of interest: 



(4.3) 


(The formula for the total cross-section differs only in that the flux is integrated only over the 
energy of bin i, rather than the whole energy range, and that we therefore do not need to divide by 
the bin width A,-.) 

We perform unfolding in these four variables using a Bayesian technique[10] with a single it¬ 
eration. The unfolding matrices Uij needed as input are predicted by our simulation. Our prediction 
for the neutrino flux d> by which we then divide is derived from a GEANT4-based simulation of 
the NuMI beamline (described further in ref. [11]). In addition, we use an in situ MINERvA mea¬ 
surement based on elastic scattering of neutrinos from atomic electrons [13] to provide a data-based 
constraint for the flux estimate. 

The cross-section vs. electron angle obtained from this procedure is given in fig. 3. We note 
that the simulation predicts substantially more events in the most forward bins than we observe in 
our data. In addition, the measured cross-section vs. shown in fig. 4a, appears to exhibit 

a noticeable migration from low to high as compared to the prediction. This is in contrast to 
the analogous cross-section measured in MINERvA using muon neutrinos, fig. 4b, which agrees 
much better—if not perfectly—with GENIE’s model. 

5. Conclusions 

Though the Vg cross-section is vitally important for neutrino oscillation searches, experimental 
challenges have prevented extensive measurement of this quantity until recently. In this preliminary 
result from MINERvA, we observe a discrepancy at low angles between the model in GENIE 2.6.2 
and our data in da/ddg. Eurthermore, we find that the Qq^- spectrum we observe appears to be 
harder for Vg CCQE than it is for v^ CCQE, in contrast to the prediction of GENIE. If substantiated 
by further study, these observations will necessitate modifications to the models currently in use in 
neutrino generators so as to ensure they correctly simulate the electron neutrino kinematics. Work 
is still ongoing to characterize the backgrounds in the Eg < 1 GeV region, and a full result will be 
published soon. 
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Figure 3: Differential cross-section vs. electron angle. Inner errors are statistical; outer are statisti¬ 
cal added in quadrature with systematic. 
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Figure 4: Differential cross-section vs. Qq^ (defined in the text) for electron and muon neutrinos. 
Inner errors are statistical; outer are statistical added in quadrature with systematic. 
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